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SANDPILES ON THE HEPTAGONAL TILING 
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Abstract. We study perturbations of the maximal stable state in a sandpile 
model on the set of faces of the heptagonal tiling on the hyperbolic plane. An 
explicit description for relaxations of such states is given. 

Dedicated to Sergei Vassilyevich Duzhin, who was a great lector and teacher - that is 
why this article is about sandpiles, one of his last courses, and passionate experimenter, 
who loved concrete calculations - that is why this article is about experimental results. 

1. Introduction 



Figure 1. The triangular tiling (black edges) and the dual hep¬ 
tagonal tiling (white edges). The hyperbolic plane is represented 
by the Klein disk model. 

Consider a Mobius triangle A C whose angles are all equal Note that such 
a triangle is unique up to an isometry of Let G be the subgroup of Aut{lP) 
generated by three reflections with respect to the sides of A. Action of G on A 


Key words and phrases, self-organized criticality, sandpiles, hyperbolic tilling. 

Research is supported in part the grant 159240 of the Swiss National Science Foundation 
as well as by the National Center of Competence in Research SwissMAP of the Swiss National 
Science Foundation. 


1 




2 


N. KALININ, M. SHKOLNIKOV 


produces the regular triangular tiling of (see [12] and Figure [^. Consider a 
graph r whose set of vertices (edges) consists of translations by G of the vertices 
(edges) of A. Note that F is a regular graph, all its vertices have degree 7. 

The dual graph of F can be embedded to by sending its vertices to the centers 
of the corresponding triangles and connecting the centers of adjacent triangles with 
a geodesic. This procedure defines the regular heptagonal tiling. We think of 
this tiling as the closest hyperbolic analog of the hexagonal honeycomb tiling of 
Euclidean plane. 

Choose a vertex O of the graph F. For any vertex of F denote by L{v) the distance 
from z) to O in F, i.e. L(v) is the minimal length of a chain of edges connecting O 
and V. For m G N denote by F^ the set of all vertices z; in F such that L{v) < m. 

For any v G consider a toppling operator T„ on the set of all integer-valued 
functions on i.e. 

Ty(j){u) = 4>{u) + A(5„(zz), 

for (j)'- Cm —t Z and u G F^. Here we denote by n{v) the set of all 7 neighbors of 
Z! in F and Sy is a function on F sending z; to 1 and vanishing at all other vertices. 
The discrete Laplacian operator A is given by 

^f{v) = -7f{v)+ 
where / is a function on F. 

A state of the sandpile model on F^ is a non-negative integer-valued function 
on Tm- We interpret a state as a distribution of sand grains (or chips) on F^. A 
toppling Ty is called legal for a state (j) if Ty(j) is a state, in other words (^(z;) > 7. 
A state (j) is called stable if there are no legal topplings for (j), i.e. (j){v) < 7 for any 
^ G T'yyi. 

Consider a state z/jq- A sequence of states z/jqi ■ • ■ i V’Af is called a relaxation of 
ipo if tpN is a stable state and for z = 0,..., A — 1 the state tpi+i is the result of 
applying a legal toppling to tpi. It is a basic statement (see um in general theory 
of abelian sandpiles that for any state z/jq there exist a relaxation and the resulting 
stable state ipN depends only on ipQ and doesn’t depend on a particular choice of 
a relaxation. Thus, we denote by ipo the resulting state of any relaxation sequence 
of ipQ. 

Denote by a function on F^ counting the number of topplings at z; in a 
relaxation of z/jq. This function is called the toppling function (or odometer) of tfo 
and for a given relaxation sequence z/jqi ■ • • i V’zv can be expressed as 

N 

'TijlO — ^ ' ^Vi ) 
i=l 

where ifi+i = Ty^ipi and Vi G Fm- The toppling function doesn’t depend on a choice 
of a relaxation sequence (see 13 El). Clearly, the result of the relaxation for ijj can 
be expressed in terms of toppling function as z/i° = z/; -|- AT^jj. 

The first version of a sandpile model was introduced in [2| where a piece of 
a standard square lattice was used instead of Fm- The model was generalized to 
arbitrary graphs in |S], the generalization essentially coincides with the chip-firing 
game, well known in combinatorics (see [13jl. The sandpile model was extensively 
studied in the case of embedded graphs including the classical square lattice and 
other tilings of Euclidean plane (see for example [D E] El 111 131313113 IHl [HI US)- 
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To our knowledge, sandpiles on hyperbolic graphs has never been described in the 
literature. Our goal is to find the analogues of our recent results (see [7]) about 
sandpiles on Euclidean lattices in the case of a hyperbolic tilling. The results 
established in this paper demonstrate that the study of sandpiles on hyperbolic 
graphs in general might be easier than on Euclidean ones. In Section]^ we discuss 
possible generalizations. 


2. Statement oe the results 

Definition 1. The maximal stable state 4>m on T^ is the state given by (j)miv) = 6 
for all V G . For a non-empty subset P of T^ consider a perturbation of the 
maximal stable state given by 

veP 

We will describe the result of the relaxation (3^ = of (j)^ and the toppling 

function . Consider a special case of the problem when P = {O}. It is easy 
visible that the states /3^ (see Figure and Pq (left side of the|^ coincide on r 2 . 
We can also observe self-similar branches on the right side of Figure Though it 



Figure 2. The state for m = 2 and P = {O}. The black cell 
is dual to the vertex O. The colors on the dual tiling represent 
numbers of sand at each vertex of r 2 . The black cell represents 0, 
dark and light grey cells represent 3 and 5 sand grains respectively. 
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Figure 3. Two pictures of the state : in the standard Klein 
model and after a homothety with ratio 0.005. 


is hard to express this fractality symbolically, we will give an easy combinatorial 
recipe for constructing the state without performing the actual relaxation (this 
is done in Proposition and Definition |^. 

Definition 2. For a state 4> denote by Mp the mass of the state given by 

M(l)= 

The quantity Mp — M(j>° represents the total number of sand grains that leave 
the system during the relaxation of p. It appears that the state looses a sub¬ 
stantial amount of sand after its relaxation and this amount doesn’t depend on the 
configuration of points P. 

Proposition 1. For all m > 0 there exist a constant Cm > 0 such that 

Cm = MpP - M/3^ 

is the same for all non-empty P C F^. Moreover, 

Cm 

hm -j-r 

m^oo [Fml 

See page for the proof. In particular, we see that in average a element of F^ 
looses more than two grains during the relaxation even if |P| = 1. Thus, adding to 
P more points produces no unstable cites and does not change the picture outside of 
the added points. This property conceptually explains the existence of a universal 
“solution” producing all the states 

Recall that L{v) denotes the distance from u to O in F. We define L{P) = 
minpgp L{p). 

Proposition 2. There exist a family Oq, oi, 02 ,... of integer-valued functions on 
the set of vertices of F with the following properties 

• Oo ^ud as coincide outside F^; 

• Og is equal to 6 on rs_i; 
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• for any m gN, non-empty P C and s = L{P) 

/3^(v) =as{v), for v G Tm\P, 

^m(p) =asip) + 1, for pGP. 

Proposition is a consequence of the following 
Theorem 1. For any non-empty P C F^ 

7^^ {v) = min(TO -I- 1 — L{v), m + 1 — L{P)). 

In other words, the toppling function essentially doesn’t depend on the number 
and the position of points P and depends only on the distance from O to P. Theorem 
and Proposition are proved on the page 

3. Combinatorics of F 

The set of vertices of F is a disjoint union of the sets F^ = Fm\Fm-i, m G Z>o. 
We call F^ the m-th level of F. For any pair of integers (m, k) we consider the set 
F^ given by 

= Hv) n f;,_i \ = k}. 

It is easy to see that Fq = Fg = {O} and F^ = F^ U F^ for m > 0. 



Figure 4. Vertices of first and second types. A type of a vertex 
is written in the circle, if it can be determined using the picture. 


We distinguish three types of vertices. The zeroth type consist of unique vertex 
O. The first type consists of vertices connected by a single edge with a previous 
level. And a vertex of the second type is connected with the previous level by two 
edges (see Figure]^. 

We are going to compute the numbers Om = |F(,^| and bm = |F^| of vertices of 
the first and second types at m-th level. It is clear that bi = 0, oi = 7, 

km+i —T and 

Hm+l T — ^m+1 T Urn- 

Therefore, the sequence 61 , 01 , 62 , 02 , 113 , 03 ,... coincides with the Fibonacci se¬ 
quence multiplied by 7. We have proven the following 

Lemma 1. For any m > 0 

6m — 702m—2 and Om — 7'02m—1; 
where Wg = 0 , Wi = 1 and Ura +2 = Hm+i + Rm- 
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Now we can compute the number of elements in Tm- 
Lemma 2. IL^I = 7u2m+i - 6. 

Proof. Indeed, the number Ir^l is equal to 

m mm 

1 + + bi) = 1 + = 1 + — Oi) = 1 + Qm+l — 0-1- 

i—1 i—1 2=1 

□ 


4. Wave action 

Consider a vertex p S r^. Let (/) be a state on r^- The “wave action” Wp(j) at p 
on (j) is defined as follows. 

Definition 3. If (j>{p) = 6 and there exists v G T^ such that v is adjacent to p and 
4>{v) = 6 then Wpcf = {TyTp(j})°. Otherwise, we define Wp4> to be (j). 

It is easy to show that the wave action doesn’t depend on the choice of v adjacent 
to p. Clearly, the states Wyip and Wp4> coincide. Therefore, for a given chain of 
vertices pi, ■ ■ ■ ,pn G such that (l>{pi) = 6 and pi+i G n{pi) all the states Wp^cj) 
coincide. 

One can use a sequence of waves to decompose the relaxation of the state (p + Sp. 
Indeed, one can verify that {(j) + dp)° = {Wp(j){v) +(5p)°. Therefore, for a sufficiently 
large N either W^(l){p) = Wp{Wp{.. .p)... ){p) < 6 and 

^^ + 5pr = W^P{p) + 5p 

or Wp (j){p) = 6 and Wp (j){v) < 6 for any v G n{p) C Tm and 

{P + SpY =Tp{W^cP{p) + 5p). 

For more details about waves see [7]. 

Now we want to describe the wave action at O applied to the maximal stable 
state (pm- Since every point v G T^ can be connected with O by a chain of adjacent 
cites and all the cites are maximal for pm, the state Wopm is the result of applying 
a toppling to pm at every point of Tm- By counting the incoming and outgoing 
grains at each vertex we conclude that Wopm is equal to 6 on Tm-i, to 2 on T^ 
and to 3 on T^ (see Figure [^. 

The main observation of this paper is that the result of applying of the second 
wave coincides with the result of applying a single wave on a smaller domain, i.e. 

WoPmiv) = Wopm-iiv) for any u G F^.i. 

Indeed, applying a toppling to Wopm at every vertex of Fm-i we get a stable state 
which is equal to 6 on Tm- 2 , to 2 on to 3 on Tm-i U and to 5 on F^ 

(see Figure]^. 

This motivates the definition of the universal family as- 

Definition 4. If s > 0 then ag is equal to 

• 6 on Fs-i, 

• 5 on Tm for m> s, 

• 3 on F^ and F^ for m > s and 

• 2 on Tl. 

The function oq is defined to be 5 on F^ and 3 on F^ for to > 0 and —1 on Fq. 


SANDPILES ON THE HEPTAGONAL TILING 


7 



Figure 5. A plot of the state Wo{4’m) (see 
the boundary of Fs. 


Definitions TjS ) 


near 



Figure 6. A plot of the state WQ(j)m near the same piece of the 
boundary as on Figure 


With this definition and the complete description of the wave action we deduce 
the following 
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Lemma 3. Consider a point p S Fq. Then + Sp and 

T^IpI (v) = niin(TO + 1 — L{p), m + 1 — L{v)). 


Proof. We decompose the relaxation process for 4>m + Sp into a sequence of waves 
actions applied to the maximal stable state. From the structure of the wave action it 
follows that the state W^cfra coincides with am-k+i restricted to Fm. In particular, 
WQ(j)m is maximal on r^-fc and is equal to 2 or 3 on = r^_^_^^\F^_^. 

Indeed, applying a toppling to WQ(f)m at each point of Tm-k we get a stable state 
which coincides with 

A point p is maximal for k <m — L{p). This implies that ITq = 

Wp^^cfrn- Using the explicit description of the wave action, we see that p is not 

maximal for and aL(p) + Sp coincides with on Combining 

these observations, we conclude that the toppling function for am + Sp is equal to 
m + 1 - L{p) on F^(p) and m - L{p) - k on F* □ 

Proof of Theorem\^ and Proposition^ Consider a configuration of vertices P C 
Fm- Let po S P be a minimizer of 


I = minL(p), 

P&p 


i.e. L{po) = 1. By Lemmathe state {(j)m + Spo)° equals to 6 exactly on Fi_i and 
therefore the state 

i'Pm+Spg)°+ ^ Sp 

pep.p^po 


is stable and is equal to f3m 


□ 


In particular, we see that some amount of sand is lost only after sending a first 
wave. We can compute the exact loss using Lemma A vertex of the first type 
on the boundary looses 4 grains of sand and a vertex of the second type looses 3 
grains. Therefore, the total loss is equal to Aom + 3bm = 7(4u2m-i + 3u2m-2) = 
7(2'U2m + ri2m+l) — 7 {pi2m + ^2m+2)- 


Proof of Proposition^ It is well known that 


Um ^ where p = i(l + v^). 


Then applying Lemma 

7p2™(l+p2) 


|r. 


Jp2m+1 


= P ^ + P 


= 75. 


□ 


5. Discussion 

The problem that we studied in this paper (sandpile on a “hyperbolic lattice”) 
is similar to the one discussed in [7] (sandpile on the two dimensional Euclidean 
lattice). For Euclidean case, the basic observation is that the amount of lost sand is 
small and therefore the resulting state is close to the maximal stable state. In fact, 
the Euclidean analogs of the states /3m (see Definition coincide with the maximal 
stable state almost everywhere. In the scaling limit, the locus of vertices where such 
a state deviates from the maximum, looks like a thin graph (or more precisely, a 
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tropical curve) passing through the perturbation points. In particular, the result 
of the relaxation strongly depends on the position of points P. Surprisingly, the 
situation is very different in the hyperbolic case. 

The key ingredient in [7] is the description of the discrete harmonic functions with 
sublinear growth on large subsets of the lattice . Also, instead of considering only 
sets like r„ we could consider polygonal subsets of Z^, it was important that the 
boundary of such subsets if defined by sets of zeroes of discrete harmonic functions. 
We do not know much about discrete harmonic functions on T, therefore we had 
only one type of subsets of T to consider. 


5.1. Scalings. Here we propose a possible procedure to obtain a scaling limit in 
the hyperbolic case. Consider a family of graphs r[n], n G N embedded to 
and satisfying the following properties. The graph T = r[l] is isohedral (i.e. tile- 
transitive, all the faces are congruent). The graph r[n -|- 1] is a result of a finite 
subdivision of each cell of r[n], moreover we ask this subdivision to be invariant 
with respect to the group of symmetries of r[n]. The set H(r[n]) of vertices of r[n] 
is dense in when n goes to infinity, i.e. lim„_,.oo dist(H(r[n]),p) = 0, for any 
p G 

Consider a geodesic-convex set H C and collection of points P C fi°. For 
simplicity assume that P C r[A^] for some fV > 0. Let (/>m be the maximal stable 
state on l/(r[TO]) C fl. Suppose that m > N, consider the perturbation = 
(0m + J2peP ^p)° of fho maximal stable. We would like to know how the limiting 
shape looks like. For example, let Cm be the subset of H(r[m]) C H where /?„ is 
not maximal. What is the Hausdorff limit (if it exists) of Cm in H? 


5.2. Boundaries. Here we propose a hyperbolic analog of the lattice polygons in 
1? with sides of rational slope. Let G be a subgroup of Aut(IHI^) consisting of all 
symmetries of the embedding of F to H^. A geodesic ^ C is called F-rational 
if there exist a point p G I and an element A S F such that p ^ Ap and I passes 
through Ap. 

It might be useful to restrict to the specific types of H. For example in [7] the 
case of polytope with rational edges plays an important role and in this case the 
results become simpler. Natural analogues of such sets would be finite intersections 
of hyperbolic half-planes whose boundaries are F-rational geodesics. 


5.3. Abstract graphs. Finally we would like to raise the main question of our 
paper in the case of general symmetric hyperbolic graphs. Let G be a finetely 
generated group and S' be a chosen finite set of its generators. For a subgroup H 
of G we construct a graph F with the set of vertices equal to C/H. Two cosets 
a,/3 G C/H are connected with an edge in this graph if the exist s G S such that 
/3 = sa. This graph is known as relative Cayley graph (see [8]). Scaling procedure 
and polygonal-type boundaries can be defined as in Sections |5.I[ |5.2| 

Let L: H(r) —)■ Z>o be the function measuring the distance from H G T, i.e. 
L{a) = distr(Q!,i7). Consider F^ = L“^([0,m]), a finite subset P of F^ and 0m 
be the maximal stable state on Tm- Is it possible to describe the relaxation for 
0m + J2peP tbe case of hyperbolic F (see [8])? 
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